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Abstract 

Modal stability parameters are extracted directly from 
aeroservoelastic flight test data by decomposition of ac- 
celerometer response signals into time-frequency atoms. 
Logarithmic sweeps and sinusoidal pulses are used to 
generate DAST closed loop excitation data. Novel 
wavelets constructed to extract modal damping and fre- 
quency explicitly from the data are introduced. The so- 
called Haley and Laplace wavelets are used to track time- 
varying modal damping and frequency in a matching 
pursuit algorithm. Estimation of the trend to aeroservoe- 
lastic instability is demonstrated successfully from anal- 
ysis of the DAST data. 


Nomenclature 


ARW-l 

first aeroelastic research wing 

C(s,u) 

continuous wavelet transform 

DAST 

drone for aerodynamic 
and structural testing 

/ 

frequency (Hz) 

m 

signal vector 

F 

Fourier transform 

g 

gravitational acceleration 

9 

window function 

9i 

time-frequency atom 

i 


K{7,4>) 

normalization constant for tla’s 
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ny 

nz 

R 

Re 

s 

tfa 

sps 

u 

Symbols 

7 _ 

T 

c 

<t> 


lateral acceleration 
normal acceleration 

residual of projection in matching pursuit 

real part 

scale 

time-frequency atom 
samples per second 
time shift 

space of parameters (s,£,u) 
space of parameters (s,— f,u) 
damping ratio 

modulation frequency of tla 
radian phase angle 
analyzing wavelet 
radian frequency 


Introduction 


Envelope expansion of new or modified aircraft often 
necessitates flight stability testing to verify safety 
margins to prevent against aeroservoelastic instabil- 
ity. Flight data is acquired for stability estimation 
and system identification to compare with analyt- 
ical predictions as a function of flight parameters 
(Mach, dynamic pressure, a). Any anomalies are 
regarded with care to guarantee safety of flight. Ex- 
citation systems are often essential to resolve sta- 
bility trends from noisy measurements since atmo- 
spheric turbulence is generally insufficient to deter- 
mine modal characteristics. 

Time-frequency representations of time varying sys- 
tems offer an alternative to traditional time and 
frequency domain stability tracking algorithms^. 
While damping estimates can be unreliable stability 
predictors^, damping and frequency trends are use- 
ful for noting changes in system dynamics as func- 
tions of flight condition. 
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This paper introduces novel wavelets for modal 
time-frequency and time-damping estimation di- 
rectly from the data, without use of intermedi- 
ate model identification schemes. Stability param- 
eters are extracted directly from the data by corre- 
lating properly designed wavelets with the system 
responses. Revisiting DAST^ ARW-1 (figure 1) 
(drone for aerodynamic and structural testing with 
the first aeroelastic research wing) data from flight 
testing performed in 1980, during which a catas- 
trophic closed loop aeroservoelastic instability was 
encountered, modal frequency and damping are es- 
timated as a function of time until instability. It is 
shown how wavelet processing can reliably identify 
time varying system stability trends. 

Stability Analysis using 
Time-Frequency Atoms 

Basis Functions 

Fourier analysis maps a signal onto an amplitude- 
frequency plane by decomposing the signal into a 
weighted sum of infinite sinusoids. Its usefulness de- 
pends on properties of time invariance, linearity, and 
stationarity of the signal. Errors are also introduced 
from time truncation and windowing schemes. Sig- 
nals that have time varying, transient, nonlinear, or 
nonstationary characteristics demand analysis func- 
tions which do not have the restrictions of Fourier 
methods. 

Wavelets are basis functions of finite duration which 
decompose a signal into scaled and time-shifted ver- 
sions of the analyzing function. The continuous 
wavelet transform is defined as^ 


atoms (tfa’s are introduced which are well lo- 
calized in time, and frequency is independent of 
scale. The tfa’s decompose a signal based on its 
local structure in terms of time, frequency, and scale 
independently. General families of tfa’s are gener- 
ated by scaling, translating, and modulating a win- 
dow function g(t) (see Appendix). Gabor wavelet 
window functions 


produce the following tfa’s over the parameter space 
of 7 = (s,£, u) (£ is frequency) 

g*y(t) = (7rs 2 )~^ e~ e J ^ 

The parameters (s,u,£) are independently adjusted 
based on local signal properties. 

A Haley wavelet is presented here as the double- 
sided exponentially damped sinusoid of figure 3 ex- x 
pressed as 

^(s,u,^) = cos(^) 

and generating the following tfa’s over j = (s,£, u) 
(see Appendix), 

5,(0 = 

V 5 

This wavelet is designed to independently extract 
modal damping and damped vibration frequency 
from a set of data, both as functions of time. There- 
fore, signals are mapped to the frequency-time and 
damping-time planes. Viscous damping ratio, C> is 
estimated empirically by the following expression 

107 rT 

c ~ * 


C{8,U) 



f(t)V(s,u,t)dt 


where u is the time translation, s the scale, and ^ 
the analyzing wavelet. Gabor wavelets are Gaussian- 
windowed sinusoids (figure 2) described as 

#($,u,t) = cos(J) 


Conventional wavelet analyses map signals onto the 
scale-time or frequency-time plane since scale and 
frequency are inversely proportional to each other. 
This restriction inhibits the extraction of modal pa- 
rameters from flight data. Hence, time-frequency 


where T is the sampling interval of the discrete time 
signal. Intuitively, as the scale, s, of the wavelet 
lengthens, damping ratio is reduced. For damping 
ratios typical of structural systems, parameter £ esti- 
mates u) n « LJdt or natural frequency u n as an equiv- 
alent approximation of damped frequency 


Finally, a Laplace wavelet is introduced as a right- 
sided version of the Haley wavelet, or a represen- 
tation of an impulse response (figure 4) and is ex- 
pressed as 


V(s,u,t) 



-I 


I cos(f) t — u > 0 
t — u < 0 
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Matching Pursuit 

The matching pursuit algorithm^ decomposes a sig- 
nal into a linear expansion of tfa’s that are from 
a redundant dictionary of functions. The wavelet 
bases are chosen according to desired features to 
be detected, or parameters to be extracted, from 
a signal. With a parametric dictionary of specific 
time-frequency atoms, matching pursuit will define 
an adaptive scale-frequency-time transform. The al- 
gorithm proceeds as follows (see Appendix): 

1. Choose initial tfa corresponding to largest co- 
efficient of C(s, u) decomposition 

2. Optimize over 7 = (s,£,u) by maximizing in- 
ner product of the signal with tfa dictionary 

3. Extract optimal tfa from the signal 

4. Repeat optimal extraction procedure until de- 
composition residual is below specified toler- 
ance 

Once the signal is decomposed to a desired degree, 
reconstruct a signal from the extracted dictionary el- 
ements. The energy distribution of the reconstructed 
signal in the frequency-time domain is generated us- 
ing a Wigner distribution^. 

Simulation Tests 

Results from the analysis of simulated data is pre- 
sented to test the accuracy of the frequency and 
damping values from application of matching pur- 
suits using the Haley and Laplace wavelets. Various 
types of simulated data are used in the tests. Haley 
decompositions are tested using the following signals 
of 100 sec each: 

1. One double-sided exponentially damped sinu- 
soid 

2. Three separated double-sided exponentially 
damped sinusoids 

3. One double-sided exponentially damped sinu- 
soid with 60% noise 

4. Three separated double-sided exponentially 
damped sinusoids with 60% noise 

Laplace decompositions are also tested using the 
above, replacing double-sided with single-sided sinu- 
soids. From the results of table 1, the capability of 
these tf a’s to extract the true values is very encour- 
aging. Values of estimated frequency and damping 
match the simulated ’true’ values quite well. For 


clean data, the estimated values are essentially ex- 
act. Noisy data results introduce a variation in some 
estimates, but the true value is always in the range 
of variation. 

Analysis of PAST Data 
with Wavelets 

The objective of the NASA DAST program was to 
pursue investigations in transonic and supersonic 
regimes using a series of aeroelastic research wings 
on drones^. These drones were equipped with 
a flutter suppression system to enable flight be- 
yond the open-loop stability boundary. Wingtip ac- 
celerometer response data at 500 sps was acquired 
for logarithmic chirps and sinusoidal doublets of a 
single cycle 20 Hz sine wave into the aileron control 
surface. 

This paper will present analysis of right wingtip ac- 
celerometer data in response to an aileron chirp in- 
put and a series of doublets just before the drone en- 
countered an aeroservoelastic instability^ with the 
flutter suppression system engaged. Several types of 
wavelet tf a’s are used in the analysis to process the 
short time data measurements of the transient modal 
responses. Each analysis (figures 5 through 8) pro- 
duces plots of the original signal (a. top), the resid- 
ual signal after the matching pursuit decomposition- 
extraction procedure (b), reconstruction of the orig- 
inal signal from the decomposition (c), frequency es- 
timates as a function of time (d), and damping esti- 
mates as a function of time, if applicable (e). 

Gabor, Haley, and Laplace wavelets are utilized with 
the matching pursuit algorithm to extract modal 
natural frequencies from the flight data. Figure 5 
presents an example of the information obtained 
from this procedure using a Gabor tfa . The top 
plot, figure 5a, shows the original accelerometer re- 
sponse from the chirp input to the aileron. Figure 5b 
presents the noise remaining in this signal after ex- 
tracting the tfa’s which result in the reconstructed 
signal of figure 5c. 

The bottom plot, figure 5d, presents the estimated 
values of the modal natural frequencies. Each cir- 
cle on this plot indicates a frequency for which the 
wavelet had a high correlation with the measured 
data. Circles representing the highest energy extrac- 
tion from the signal by the tfa correspond to modal 
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responses. The increasing frequency as a function of 
time for the wavelet correlations agrees with the log- 
arithmically increasing frequency in the chirp signal 
and demonstrates the additional time information 
that cannot be obtained from traditional Fourier 
analysis. 

Modal frequencies similar to figure 5d are computed 
using the Haley and Laplace wavelet tfa’s . These 
frequencies agree closely with the Gabor wavelet 
analysis; however, the Haley and Laplace tfa’s ad- 
ditionally extract modal damping. For this reason, 
the remaining discussion will consider only the Haley 
and Laplace wavelet analysis. 

The matching pursuit algorithm is able to extract 
information about multiple modal responses by cor- 
relating wavelets of different frequency and damping 
(related to wavelet scale) with the measured data. 
Modal estimates are separated for clarity of pre- 
sentation with figure 6d-e displaying the estimated 
Wing Torsion mode and figure 7d-e displaying the 
estimated Wing Bending modal estimates. 

Frequency estimates of the Wing Torsion mode in 
figure 6d remain nearly constant at 25.5 Hz through- 
out the data set. The damping value of £ « 0.075 
also shows little variation throughout the series of 
doublets. These values agree with the results of pre- 
vious post- flight estimates! 1 * 7 J. 

The Wing Bending mode is of particular interest 
since this mode went unstable and resulted in the 
instability which destroyed the drone. The infor- 
mation obtained by the matching pursuit analysis 
with the Haley wavelet is given in figure 7. The 
estimated natural frequency for this mode shows 
some variation at different test points throughout 
the data set but remains centered around 20 Hz. Es- 
timated damping, however, shows significant varia- 
tion. The tfa analysis tracks modal damping de- 
creasing steadily from C ~ 0.02 to near zero after 40 
seconds of data. 

Similary, the Laplace tfa analysis given in figure 8 
indicates the impending instability. The estimated 
damping in figure 8e shows the same decrease in 
damping from C ~ 0.02 to near zero after 40 seconds 
of data. 

The Haley and Laplace wavelet analysis agree with 
the true behaviour of the DAST drone. Each wavelet 


indicates the drone to be on the verge of instability 
for the Wing Bending mode after 40 seconds of data. 
In reality, this mode did go unstable immediately 
after the 40 seconds during which this data set was 
recorded. 

Algorithms are currently being developed which im- 
plement this matching pursuit strategy with reduced 
computational effort. The large amount of energy re- 
maining in the signals after the correlated wavelets 
are extracted in figures 6b, 7b, and 8b is due to ter- 
minating the analysis before all the energy could be 
extracted by the wavelets. The process was stopped 
once the main features of the data were determined. 

Conclusions 

An aeroservoelastic system was used to test a 
method for time-frequency estimation of modal pa- 
rameters using wavelets. Time-frequency-scale data 
processing allows improved visualization and under- 
standing of the signal information content. Tuned 
basis functions for specific applications assist in iden- 
tifyingstructure (shape, amplitude, frequency, dura- 
tion, and timing of events) in the data. 

Wavelets are exploited to parameterize the data such 
that stability trends can be detected accurately in a 
time-varying flight test scenario. In this appication, 
time-frequency atoms are used in a matching pursuit 
algorithm to extract modal frequency and damping 
as functions of time. DAST data of the final seconds 
before a closed loop instability are analyzed to vali- 
date this approach as a stability indicator and modal 
tracker. Future work will enhance this capability for 
application in an on-line test environment. 
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Time-frequency Atoms 


Time-frequency atoms are functions well localized in 
both time and frequency which can decompose sig- 
nals adaptively according to particular signal struc- 
tures. General families of time-frequency atoms 
can be generated by scaling, translating and mod- 
ulating a window function g(t) £ L 2 (R). Assume 
g(t) is real, continuously differentiable, || g(t) || = 1 , 
SToo 0(0 7 ^ 0’ an d 5^(0) 7^ 0. For any scale s > 0, fre- 
quency modulation f, and time translation u, denote 
7 = (s,(,u)er = R + x R 2 and define 

g,(t) = ±g(i^L)e^. (1) 

The factor ^ normalizes || g(t) || = 1. Generally 
g(t) is even, centered at u with most of its energy 
concentrated near u, and determined by size pro- 
portional to s. The Fourier transform of g(t) is 

g y ( u ) = ^g(s(u-0)e- i(u, ~ ()u 

For even g(t ), g(u) is even and centered at frequency 
u) = f. Its energy is concentrated near center and 
modulation frequency £. 

To analyze signal structures of varying sizes, it is 
necessary to use time-frequency atoms of different 
scales. A wavelet transform decomposes signals 
over different scales, but the resulting atoms are re- 
stricted by the relation f a j. For signals that in- 
clude scaling and highly oscillatory structures, the 
appropriate constraints on the scale and modula- 
tion parameters of the time-frequency atoms p 7 can- 
not be defined apriori. Hence, an adaptive selec- 
tion of the elements of the dictionary V = (<77 n ) 7e r 
is needed depending on the local properties of sig- 
nal f(t). Finite linear expansions of time- frequency 
atoms are dense in L 2 (R), so the dictionary V is 
complete. 

Matching Pursuit Algorithm 


where Rf is the residual vector after projecting / 
onto vector g lo . Since g lo is orthogonal to Rf , 

II /II 2 = I < f<9y 0 > | 2 + II Rf II 2 - 

Minimize || Rf || by choosing g lQ £ V such that | < 
/, g lo > | is maximum, or at least find a vector g lo 
that is sub-optimal in the sense that 

\<f,9t 0 >\ > A sup | < /, > | 

7<=r 

where 0 < A < 1 is an optimality parameter. Now, 
by induction, suppose R°f = f so at the n th step 
the residue R n f is matched 

I < Rn f> 9-t„ > I > A sup | < R n f, g y > |. 

7 eT 

Now decompose R n f into 

R n f = <R n f,9i n >9i n + R n+l f 

defining the residue at step n -f 1. Since R n+l f is 
orthogonal to g ln , 

\\ Rn f II 2 = I < R n fi9y„ > I 2 + II R n+1 f || 2 . 

and the decomposition up to the m th stage results 
in 

m — 1 

/ = J2( Rn f ~ R n+1 f) + Rm f 

n = 0 
m-1 

= X < R n f>9- 1n >9y„ + R m f- 

n=0 

Therefore, the original signal / is decomposed into a 
sum of the dictionary elements chosen to best match 
its residues. Although this decomposition is nonlin- 
ear, energy conservation is preserved since 

m— 1 

|| /|| 2 = J 2 \ < Rn f ’ 9 i n >\ 2 + \\ R m fW 2 - 

n = 0 


To compute a linear expansion of a signal / £ L 2 (R) 
from the set V , successive approximations of / are 
performed with orthogonal projections on the vec- 
tors in V to best match dominant patterns in the 
data. The following algorithm, called matching pur- 
suit, performs the adaptive decomposition. 

A signal vector / can be decomposed into (<, > de- 
notes inner product) 

/ = < f,9t 0 > 9y 0 + Rf 


Gabor Atoms 


The scale modulated form of the complex Gabor 
atom of equation 1 can be written (£ = 2 tt/ c , where 
f c is the center modulation frequency in Hertz) 


9-r(t) 

9y~(t) 


<,-*( *=*)%•'«* 


£± e -iW e -*t 

\A 
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where 7 = (s,£,u) and 7“ = (s,-£,u), respectively. 
The C 0 to get || g y (t) || = 1 is C 0 = 7r“« , so 

g y (t) = (7rs 2 )“^e _ ^~^e , ^ t . 

When a signal f(t) is real, dictionaries of real time- 
frequency atoms must be used to get a decom- 
position with real expansion coefficients. For any 
7 = (s,£,u), with f ± 0 and any phase <j) E [0, 27t), 
define 


5(7,0) = V 1 ) cos(ft + 0). 

The phase <j> now appears explicitly rather than im- 
plicitly as in the complex representation of equa- 
tion 1. It is the phase angle of the inner product 
< R n f,g y > for each term of the decomposition. In 
the real Gabor case, 


9 ( 7,« 


e a( . )’ cos(£f + <£). 
\/s 


where real atoms are related to complex atoms by 

5(7,0) = ^ftW + e-VW) 

Real matching pursuits are not equivalent to com- 
plex decompositions because <7 7 (t)and </ 7 -(<)are 
not orthogonal. The constant K^^to maintain 
|| <7( 7(< £) || = 1 for real atomic decompositions be- 
comes 


A(7>) \] 1 + Re(e‘ 2 ^ < g y (t),g y -(t) » 

and for real Gabor atoms^ 


( 2 ) 


K, 


V 1 + Ee(e ,2 (^ + ^ u )e - ( 5 G 3 ) ' 


Haley Atoms 


to -p, 

<»,((),»,-«)> = .-’IW* 1 ® 

= i-P 3 e-ai^le-'^dp 

s J “ cc 

which is also j multiplied by the Fourier transform 

of c -2| £ f i | w ith w replaced by 2£. Since T[c~ 2t ] = 
, from the Fourier scaling and time-shift prop- 
erties, 

^•[ar(a<)] 

= » [ j+ifcr 

and 

= e - iwU X(f) 

jr[ e -2|^l] _ e’“ u X(w) 

4se ,w 
4 -f (so;) 2 ' 

Multiplying this result by j and replacing w by 2f 
yields 

<5 7 (0.57-(0> = 7 

5 


1 + («0 2 

for the expression in equation 2. Hence, for the Ha- 
ley atom, 


45 e ’ 2 <» 
[ 4+4£ 3 j 3 

e .*2C« 


-X ( L ) 

a v a ' 
= X{u) 




2 




1 + Re ( 


■hW ) 


Laplace Atoms 


Laplace atoms are a trivial modification of Haley 
atoms by a factor of resulting in -f- = 1 to get 
|| g y || = 1, so that 


In the so-called Haley atomic decomposition intro- 
duced in this paper, the window function is #(2) = 
C 0 e“M, and C 0 = 1 to get || g y || = 1, so that 

j,«) = 


Sj(t) = sj^e 
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K±JI 

.006, .008, .004 

.005 


\ Laplace Wavelet Results 
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Table 1: Results of Matching Pursuit Tests using Sim- 
ulated Data - / is frequency (Hz), Q is damping ratio 



Figure 1: DAST ARYV-1 



Figure 2: Gabor wavelet 



Figure 3: Haley wavelet 


Mode 

Calculated | Measured 

Symmetric Mode 

1 1 

Wing 1" Bending 

9.1 

9.6 

Fuselage 1** Bending 

16.5 

16.2 

Wing Bending-torsion 

29.6 

29.1 

Antisymmetric Mode 

i — 

Fuselage 1 * ' Bending 

12.3 

13.5 

Wing l* 1 Bending 

21.7 

19.3 

Stabilizer 1 * 1 Bending 

30.0 

27.0 


Table 2: DAST ARW-1 Calculated and GVT-measured 
Elastic Modal Frequencies (Hz) 


- 
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Figure 4: Laplace wavelet 
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